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Amendments to the Claims: 

This listing of claims will replace all prior versions, and listings, of claims in the 

application: 

1. (Original) A method for reconstructing an image of a scattering medium, comprising: 

directing energy into the scattering medium at a source location on the 
scattering medium; 

measuring the energy emerging from the scattering medium at a detector 
location on the scattering medium; 

selecting an initial guess of internal properties of the scattering medium; 
predicting the energy emerging from the scattering medium using an 
equation of radiative transfer, wherein the prediction is a function of the 
initial guess; 

generating an objective function based on a comparison of the prediction 
with the measurement; 

generating a gradient of the objective function by a method of adjoint 
differentiation; 

modifying the initial guess of the properties of the scattering medium 
based on the gradient of the objective function; and 
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generating an image representation of the internal properties of the 
scattering medium. 

2. (Presently Amended) The method according to claim 1, further comprising repeating the 
predicting of the energy emerging from the scattering medium based on the modified initial 
guess, generating the objective function and modifying the initial guess, until at least one of a 
predetermined number of repetitions has occurred and— or the objective function reaches a 
predetermined threshold. 

3. (Presently Amended) The method according to claim 1, wherein the prediction depends 
on the boundary conditions. 

4. (Presently Amended) The method according to claim 3, wherein the boundary conditions 
account for a refractive index mismatch at an interface between the medium and at least one ef 
the detectors-and or source. 

5. (Original) The method according to claim 1, wherein the prediction comprises an 
iterative process producing intermediate results. 

6. (Original) The method according to claim 5, wherein the intermediate results of the 
prediction are stored. 

7. (Original) The method according to claim 6, wherein generating the gradient of the 
objective function by adjoint differences uses the intermediate results of the prediction. 
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8. (Original) The method according to claim 7, wherein generating the gradient comprises 
stepping backward through the intermediate results of the prediction. 

9. (Original) The method according to claim 1, wherein the equation of radiative transfer is 
time independent. 

10. (Original) The method according to claim 9, wherein the time independent equation of 
radiative transfer is: 

where ^(r,^) is the radiance at the spatial position r directed into a unit solid angle co, S(r,w) is 
the energy directed into the medium at spatial position r into a unit solid angle co, /u s is the 
scattering coefficient, /u a is the absorption coefficient and p(co,of) is the scattering phase 
function. 

1 1 . (Original) The method according to claim 1 0, wherein the scattering phase function is: 

P(COS0) = ZTT 

PK } 2(l + g 2 -2gcos0) 3/2 
where 6 is the angle between the two unit solid angles co and co 
and g is the anisotropy factor. 

12. (Original) The method according to claim 1, wherein the equation of radiative transfer is 
time dependent. 
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13. (Original) The method according to claim 12, wherein the time dependent equation of 
radiative transfer is: 



where ^(r^t) is the radiance at the spatial position r directed into 
a unit solid angle cy, S(r,w, t) is the energy directed into the 
medium at spatial position r into a unit solid angle co, fu s is the 
scattering coefficient, /u a is the absorption coefficient and p{o),co f ) 
is the scattering phase function. 



14. (Original ) The method according to claim 13, wherein the scattering phase function is: 



15. (Presently Amended) The method according to claim 1, wherein the properties include at 
least one of a scattering coefficient, an absorption coefficient, an anisotropy factor, aftd-or a 
scattering phase function. 

16. (Presently Amended) The method according to claim 1, wherein the objective function is 
a normaliz e d comparison of the predicted energy and the measured energy. 



1 ffF(r,cM) 
c dt 





where 9 is the angle between the two unit solid angles co and a> 



and g is the anisotropy factor. 
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17. (Presently Amended) The method according to claim 1, wherein the objective function is 
based on the normaliz e d sum of the differences between the predicted energy and the measured 
energy for each source detector pair, wherein a source detector pair is formed between each 
source location and each detector location. 

18. (Original) The method according to claim 1, wherein the objective function is: 

<p=-z(p i -M i y 

z i 

where M; represents the actual measurements and the Pi represents the predicted 
measurements for each source defector pair i, m is the number of source detector pairs, where a 
source detector pairs is formed between each source location and each detector location. 

19. (Original) The method according to claim 1, further comprising minimizing the objective 
function. 

20. (Original) The method according to claim 19, wherein minimizing the objective function 
includes a one dimensional line search. 

21. (Original) The method according to claim 20, wherein the one dimensional line search is 
performed along a direction of the gradient of the objective function. 

22. (Original) The method according to claim 20, wherein the one dimensional line search is 
performed along a gradient-dependent direction. 
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23. (Original) The method according to claim 1, wherein the energy comprises near infra-red 
energy. 

24. (Original) The method according to claim 1, wherein the scattering medium contains 
regions wherein the scattering coefficients are not substantially greater than the absorption 
coefficients. 

25. (Original) The method according to claim 1, wherein the scattering medium contains a 
low scattering region embedded in a high scattering region. 

26. (Original) The method according to claim 1, wherein the predicted energy is determined 
using finite element methods. 

27. (Original) The method according to claim 1, wherein the predicted energy is determined 
using finite difference methods. 

28. (Presently Amended) A method for imaging the spatial optical properties of tissue, 
comprising: 

(a) directing energy into the scattering medium at a source 
location on the tissue; 

(b) measuring the energy emerging from the scattering medium 
at a detector location on the tissue; 
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(c) selecting and an initial guess of the spatial optical 
properties of the tissue; 

(d) predicting the energy emerging from the tissue using an 
equation of radiative transfer in an iterative process, wherein the 
prediction is a function of the initial guess and a refraction index 
mismatch at a boundary of the tissue, and the iterative process 
generates a plurality of intermediate predictions; 

(e) generating an objective function based on a 
normaliz e d comparison of the prediction with the measured energy 
emerging from the scattering medium; 

(f) generating a gradient of the objective function by adjoint 
differentiation; 

(g) modifying the initial guess of the spatial properties of the 
tissue based on the gradient of the objective function; 

(h) repeating steps (d) through (g) until at least one of a 
threshold of modifications to the initial guess is reached and-or the 
objective function reaches a threshold; and 

(j) generating an image representation of the spatial optical 
properties of the tissue. 
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29. (Presently Amended) A system for reconstructing an image of a scattering medium, 
comprising: 

a source for directing energy into the scattering medium at a^source location on the 
scattering medium; 

a detector for measuring the energy emerging from the scattering medium at a detector 
location on the scattering medium; 

an initial guess of internal properties of the scattering medium; 

means for predicting the energy emerging from the scattering medium using an equation 
of radiative transfer, wherein the prediction is a function of the initial guess; 

means for generating an objective function based on a comparison of the prediction with 
the measurement; 

means for generating a gradient of the objective function by a method of adjoint 
differentiation; 

means for modifying the initial guess of the properties of the scattering medium based on 
the gradient of the objective function; and 

means for generating an image representation of the internal properties of the scattering 
medium. 
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30. (Presently Amended) The system according to claim ±29, further comprising means for 
repeating the predicting of the energy emerging from the scattering medium based on the 
modified initial guess, generating the objective function and modifying the initial guess, until at 
least one of a predetermined number of repetitions has occurred and-or the objective function 
reaches a predetermined threshold. 

31. (Presently Amended) The system according to claim 4-29, wherein the prediction 
depends on the-boundary conditions. 

32. (Presently Amended) The system according to claim 31, wherein the boundary 
conditions account for a refractive mismatch at an interface between the medium and at least one 
effee-detectors-and-or source. 

33. (Presently Amended) The system according to claim ±29, wherein the prediction 
comprises an iterative process producing intermediate results. 

34. (Original) The system according to claim 33, wherein the intermediate results of the 
prediction are stored. 

35. (Presently Amended) The system according to claim 34, wherein generating the gradient 
of the objective function by adjoint diff e r e nc e s differentiation uses the intermediate results of the 
prediction. 
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36. (Original) The system according to claim 35, wherein generating the gradient comprises 
stepping backward through the intermediate results of the prediction. 

37. (Presently Amended) The system according to claim ±29, wherein the equation of 
radiative transfer is time independent. 

38. (Original) The system according to claim 37, wherein the time independent equation of 
radiative transfer is: 

where ^(r,^) is the radiance at the spatial position r directed into a unit solid angle to, S(r,w) is 
the energy directed into the medium at spatial position r into a unit solid angle at, /u s is the 
scattering coefficient, ju a is the absorption coefficient and p{o),co r ) is the scattering phase 
function. 

39. (Original) The system according to claim 38, wherein the scattering phase function is: 

^ 2(l + g 2 -2gcos0) 3/2 

where 0 is the angle between the two unit solid angles co and 0) and g is the 
anisotropy factor. 

40. (Presently Amended) The system according to claim ±29, wherein the equation of 
radiative transfer is time dependent. 

11 

912224 v2 



Patent 

Serial No. 09/767,230 
Docket No. 0887-4150US1 



41. (Original) The system according to claim 40, wherein the time dependent equation of 
radiative transfer is: 



where ^(r^t) is the radiance at the spatial position r directed into a unit solid angle co, S(r,w, t) 
is the energy directed into the medium at spatial position r into a unit solid angle co, ju s is the 
scattering coefficient, ju a is the absorption coefficient and pico.co") is the scattering phase 
function. 

42. (Original) The system according to claim 41 , wherein the scattering phase function is: 



where 6 is the angle between the two unit solid angles co and co ', and g is the 
anisotropy factor. 

43. (Presently Amended) The system according to claim ±29, wherein the properties include 
at least one of a scattering coefficient, an absorption coefficient, an anisotropy factor, and-or a 
scattering phase function. 

44. (Presently Amended) The system according to claim 429, wherein the objective function 
is a normaliz e d comparison of the predicted energy and the measured energy. 



1 &¥(r,(Q,t) 
c dt 




p(cos0) = 




2(l + g 2 -2gcos6>) 



3/2 
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45. (Presently Amended) The system according to claim 4-29, wherein the objective function 
is based on the normaliz e d sum of the differences between the predicted energy and the measured 
energy for each source detector pair, wherein a source detector pair is formed between each 
source location and each detector location. 

46. (Presently Amended) The system according to claim 429, wherein the objective function 
is: 

1 m , o 

z i 

where Mi represents the actual measurements and Pi represents the predicted 
measurements for each source detector pair, m is the number of source detector pairs, where a 
source detector pairs is formed between each source location and each detector location. 

47. (Presently Amended) The system according to claim 4-29, further comprising minimizing 
the objective function. 

48. (Original) The system according to claim 47, wherein minimizing the objective function 
includes a one dimensional line search. 

49. (Original) The system according to claim 48, wherein the one dimensional line search is 
performed along a direction of the gradient of the objective function. 
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50. (Original) The system according to claim 49, wherein the one dimensional line search is 
performed along a gradient-dependent direction. 

51. (Original) The system according to claim 50, wherein the energy comprises near infra- 
red energy. 

52. (Presently Amended) The system according to claim 4-29, wherein the scattering medium 
contains regions wherein the scattering coefficients are not substantially greater than the 
absorption coefficients. 

53. (Presently Amended) The system according to claim ±29, wherein the scattering medium 
contains a low scattering region embedded in a high scattering region. 

54. (Presently Amended) The system according to claim 4-29, wherein the predicted energy 
is determined using finite element methods. 

55. (Presently Amended) The system according to claim ±29, wherein the predicted energy 
is determined using finite difference methods. 

56. (Presently Amended) A system for imaging the spatial distribution of optical properties 
of tissue, comprising: 

(a) a source for directing energy into the scattering medium at a source location on 
the tissue; 
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(b) a detector for measuring the energy emerging from the scattering medium at a 
detector location on the tissue; 

(c) an initial guess of spatial optical properties of the tissue; 

(d) means for predicting the energy emerging from the tissue using an equation of 
radiative transfer in an iterative process, wherein the prediction is a function of the initial guess 
and a refraction index mismatch at a boundary of the tissue, and the iterative process generates a 
plurality of intermediate predictions; 

(e) means for generating an objective function based on a normaliz e d comparison of 
the prediction with the measured energy emerging from the scattering medium; 

(f) means for generating a gradient of the objective function by adjoint 
differentiation; 

(g) means for modifying the initial guess of the spatial properties of the tissue based 
on the gradient of the objective function; 

(h) means for repeating steps (d) through (g) until at least one of a threshold of 
modifications to the initial guess is reached and-or the objective function reaches a threshold; and 

(j) means for generating an image representation of the spatial optical properties of 
the tissue. 

57. (Original) Computer executable software code stored on a computer readable medium, 
the code for reconstructing an image of a scattering medium, comprising: 
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code to direct energy into the scattering medium at a source 
location on the scattering medium; 

code to measure the energy emerging from the scattering medium 
at a detector location on the scattering medium; 
code to receive an initial guess of internal properties of the 
scattering medium; 

code to predict the energy emerging from the scattering medium 
using an equation of radiative transfer, wherein the prediction is a 
function of the initial guess; 

code to generate an objective function based on a comparison of 
the prediction with the measurement; 

code to generate a gradient of the objective function by a method 
of adjoint differentiation; 

code to modify the initial guess of the properties of the scattering 
medium based on the gradient of the objective function; and 
code to generate an image representation of the internal properties 
of the scattering medium. 
58. (Presently Amended) Computer executable software code stored on a computer readable 
medium, the code for imaging the spatial distribution of optical properties of tissue, comprising: 

(a) code to direct energy into the scattering medium at a source location on the tissue; 
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(b) code to measure the energy emerging from the scattering medium at a detector 
location on the tissue; 

(c) code to receive an initial guess of spatial optical properties of the tissue; 

(d) code to predict the energy emerging from the tissue using an equation of radiative 
transfer in an iterative process, wherein the prediction is a function of the initial guess and a 
refraction index mismatch at a boundary of the tissue, and the iterative process generates a 
plurality of intermediate predictions; 

(e) code to generate an objective function based on a normaliz e d comparison of the 
prediction with the measured energy emerging from the scattering medium; 

(f) code to generate a gradient of the objective function by adjoint differentiation; 

(g) code to modify the initial guess of the spatial properties of the tissue based on the 
gradient of the objective function; 

(h) code to repeat steps (d) through (g) until at least one of a threshold of modifications to 
the initial guess is reached and-or the objective function reaches a threshold; and 

(j) code to generate an image representation of the spatial optical properties of the tissue. 
59. (Original) A computer readable medium having computer executable software code 
stored thereon, the code for reconstructing an image of a scattering medium, comprising: 

code to direct energy into the scattering medium at a source 
location on the scattering medium; 
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code to measure the energy emerging from the scattering medium 
at a detector location on the scattering medium; 
code to receive an initial guess of internal properties of the 
scattering medium; 

code to predict the energy emerging from the scattering medium 
using an equation of radiative transfer, wherein the prediction is a 
function of the initial guess; 

code to generate an objective function based on a comparison of 
the prediction with the measurement; 

code to generate a gradient of the objective function by a method 
of adjoint differentiation; 

code to modify the initial guess of the properties of the scattering 
medium based on the gradient of the objective function; and 
code to generate an image representation of the internal properties 
of the scattering medium. 
60. (Presently Amended) A computer readable medium having computer executable 
software code stored thereon, the code for imaging the spatial distribution of optical properties of 
tissue, comprising: 

(a) code to direct energy into the scattering medium at a source location on the tissue; 



912224 v2 



18 



Patent 

Serial No. 09/767,230 
Docket No. 0887-4150US1 

(b) code to measure the energy emerging from the scattering medium at a detector 
location on the tissue; 

(c) code to receive an initial guess of spatial optical properties of the tissue; 

(d) code to predict the energy emerging from the tissue using an equation of radiative 
transfer in an iterative process, wherein the prediction is a function of the initial guess and a 
refraction index mismatch at a boundary of the tissue, and the iterative process generates a 
plurality of intermediate predictions; 

(e) code to generate an objective function based on a normaliz e d comparison of the 
prediction with the measured energy emerging from the scattering medium; 

(f) code to generate a gradient of the objective function by adjoint differentiation; 

(g) code to modify the initial guess of the spatial properties of the tissue based on the 
gradient of the objective function; 

(h) code to repeat steps (d) through (g) until at least one of a threshold of modifications to 
the initial guess is reached and-or the objective function reaches a threshold; and 

(j) code to generate an image representation of the spatial optical properties of the tissue. 
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